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ABSTRACT 

We derive the evolution of the energy deposition in the intergalactic medium (IGM) 
by dark matter (DM) decays/annihilations for both sterile neutrinos and light dark 
matter (LDM) particles. At 2; > 200 sterile neutrinos transfer a fraction /abs ^ 0.5 of 
their rest mass energy into the IGM; at lower redshifts this fraction becomes <, 0.3 
depending on the particle mass. The LDM particles can decay or annihilate. In both 
cases /abs ~ 0.4 — 0.9 at high (> 300) redshift, dropping to « 0.1 below z — 100. 
These results indicate that the impact of DM decays / annihilations on the IGM thermal 
and ionization history is less important than previously thought. We find that sterile 
neutrinos (LDM) decays are able to increase the IGM temperature by z = 5 at most 
up to 4 K (100 K), about 50-200 times less than predicted by estimates based on the 
assumption of complete energy transfer to the gas. 

Key viTords: cosmology: dark matter - neutrinos - ISM: evolution 



1 INTRODUCTION 



According to 3-yr WMAP ( Wilkinson Microwave Anisotropy 
Probe) results (Spergel et al. 2006), the dark matter (DM) 
constitutes about 20% of the cosmic energy density. How- 
ever, the nature of such elusive component remains unclear. 

As proposed DM candidates might induce drastically 
different evolutionary effects depending on their properties 
(e.g. velocity dispersion) , one hopes to be able to select suit- 
able DM candidates from the comparison between their pre- 
dicted effects and observations. 

Potentially important cosmological effects might be in- 
duced if DM particles either decay or annihilate, as predicted 
by fundamental physics theories. For example, sufficiently 
light DM particles (mass < 100 MeV; Boehm et al. 2004; 
Ascasibar et al. 2006) can annihilate, or decay into lighter 
particles (Hooper & Wang 2004; Picciotto & Pospelov 2005; 
Ascasibar et al. 2006) remaining good DM candidates. The 
products of DM decays/annihilations can be photons, neu- 
trinos, electron-positron pairs, and/or more massive parti- 
cles, depending on the mass of the progenitor. 

The decay /annihilation of DM particles into e"*" — e~ 
pairs has been recently invoked to explain the observation, 
by the SPI spectrometer aboard ESA's INTEGRAL sateb 
lite, of an excess in the 511-keV line emission from the Galac- 
tic bulge (Knodlseder et al. 2005). Although exotic, this idea 
has triggered many theoretical studies (Boehm et al. 2004; 
Hooper & Wang 2004; Picciotto & Pospelov 2005; Kawasaki 
& Yanagida 2005; Kasuya & Takahashi 2005; Casse & Fayet 



2006; Ascasibar et al. 2006), aimed at constraining DM prop- 
erties through SPI/INTEGRAL observations. 

The products of decays and/or annihilations are ex- 
pected to interact with the intergalactic medium (IGM), 
transferring part of their energy. If so, DM decays and 
annihilations might change the IGM thermal/ionization 
history in a sensible and detectable way. Various flavors 
of this mechanism have been investigated in a consider- 
able number of studies (Hansen & Haiman 2004; Chen 
& Kamionkowski 2004; Pierpaoli 2004; Padmanabhan & 
Finkbeiner 2005; Mapelli & Ferrara 2005; Biermann & 
Kusenko 2006; Mapelli, Ferrara & Pierpaoli 2006, hereafter 
MFP06; Zhang et al. 2006). 

The above studies, though, made the simplifying as- 
sumptions that either (i) the energy injected by DM de- 
cays/annihilation is entirely absorbed by the IGM (Hansen 
& Haiman 2004; PierpaoU 2004; Biermann & Kusenko 2006; 
MFP06), or (ii) leave the absorption efficiency as a free 
parameter (Padmanabhan & Finkbeiner 2005; Zhang et 
al. 2006), or (iii) make a partial treatment of the physical 
processes responsible for the energy redistribution (Chen & 
Kamionkowski 2004; MapeUi & Ferrara 2005). 

In this paper we model in detail, for the first time, 
the physical processes governing the interaction between the 
IGM and the decay/annihilation products, and we derive the 
fraction of energy actually absorbed (Section 2). We restrict 
our analysis to the case in which the decay /annihilation 
products are photons, electron-positron pairs, or neutrinos 
(which are assumed to have negligible interactions with mat- 
ter), because of the uncertainties in modeling the cascade 
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associated with more massive product particles. For pho- 
tons (Section 2.1) we include the effects of Compton scat- 
tering and photo-ionization; for pairs, the relevant processes 
are inverse Compton scattering, collisional ionizations, and 
positron annihilations (Section 2.2). 

Our model exhaustively describes the behavior of rel- 
atively light DM candidates (i.e. warm and cold DM par- 
ticles with mass lower than ~ 100 MeV), whose only de- 
cay/annihilation products are photons, pairs and neutrinos. 
As an application, wo explore the effects of sterile neutrinos 
(mass 2-50 keV) and of viable light dark matter particles 
(LDM; mass 1-10 MeV) on the IGM thermal and ionization 
history (Section 3). 

We adopt the best-fit cosmological parameters after the 
3-yr WMAP results (Spergel et al. 2006), i.e. Vl^ = 0.042, 

= 0.24, Odm = - t^b = 0.198, = 0.76, h = 0.73, 
Ho = 100 /i km s"^ Mpc-^ 



2 ENERGY INJECTION IN THE IGM 

We consider the case of DM decays and annihilations result- 
ing in the direct production of photons, electron-positron 
pairs, and/or neutrinos, and study the absorption of the en- 
ergy of such particle products by the IGM. 

Some fraction of the energy of the newly created parti- 
cles is immediately absorbed; the remainder will remain in 
the form of a background, eventually absorbed at later times. 

In practice, at each redshift z, it is convenient to dis- 
tinguish the particles which wore produced by DM decays 
and/or annihilations according to their "production red- 
shift" z' and energy E. We define n{z, z') as the number 
(per baryon) of product particles which at redshift z can 
still inject energy in the IGM, and were produced at red- 
shift z' > z. Their energy spectrum is ^{z, z' , E). 

The rate of energy absorption per baxyon in the IGM 
at redshift z is the sum of the contributions of all "prod- 
uct particles", obtained by the integration over production 
redshift and energy 



e{z) 



dz' j dE^{z,z',E)E <t>iz, E) 



(1) 



where cl>{z, E) is the fraction of the energy i? of a parti- 
cle which is absorbed by the IGM per unit time, calcu- 
lated at redshift z. We assume 0max=llOO, which is approx- 
imately the redshift of the last scattering surface (Spergel et 
al. 2006). 

However, using equation (1) is both complicated and 
unnecessary: in fact, we know that the energy spectrum of 
"fresh" DM decay /annihilation products is essentially mono- 
energetic, and it is reasonable to expect that it will remain 
peaked at the average energy. Therefore, we assume that the 
energy integration inside equation (1) can be safely elimi- 
nated by using the average energy E[z, z') of the particles 



SO that equation (1) becomes 
^(^) = / dz n{z,z) E{z,z')(j){z,E{z,z')). 

J z 

The evolution of n{z,z'), E{z,z'), and (j){z,E) obviously 



(2) 



(3) 



depends on the type of DM decay /annihilation product 
we are considering. If the product particles are neutrinos 
4>,^{z,E) = 0; the cases of photons and electron-positron 
pairs are discussed in the following subsections. 



2.1 Photons 

If the DM decays or annihilations result in the pro- 
duction of photons, the only two important energy loss 
mechanisms in the considered energy and redshift ranges 
(25 eV < S < 50 MeV; z < 1000 ) are ionizations and Comp- 
ton scattering on cold matter (see Zdziarski & Svensson 
1989; hereafter ZS89). All the energy lost by the photons 
is absorbed by the IGM, so we can write 



.{z,E) + < 



.{z,E). 



(4) 



The photo-ionization term can be expressed as 

J, I TP\ 0-Hc+H(i?) , ^ 

0T,,ion(^, E) = — Nh{z) C, 

where Nb{z) ~ 2.5 x lO^^l + zfi^lbh^ /0.0224) cm"^ is the 
number density of the baryons at redshift z, and 



(5) 



auc+niE) = CTHc + 12crH ~ 5.1 x 10 



V250 eV/ 



cm (6) 



(with p = 3.30 for E > 250 eV, p = 2.65 for 25eV < < 
250 eV) is the photoionization absorption cross-section per 
helium atom of the cosmological mixture of H and He (see 
equation 3.2 of ZS89). Equation (5) implicitly assumes that 
the IGM is mostly neutral, which is true between the hy- 
drogen recombination at 2; ~ 1100, and its reionization at 
« ~ 6 — 15. 

The Compton scattering term is 



,{z,E) = [cTT^giO]N4z)c 



(7) 



where ^ = E / {meC?), Na{z) is the electron number density, 
and the product (Tt C ff(C) gives the average fraction of en- 
ergy which is lost by a photon for each electron on its path 
(ctt — 6.65 X 10~^^ cm~^ is the Thomson cross-section; see 
equation 4.9 of ZS89 for the definition of the function 5(^)). 
In equation (7) the electron number density accounts for 
both free and bound electrons, as the energy losses due to 
Compton scattering only become important when E is so 
high that the interaction is insensitive to whether an elec- 
tron is bound or free (Chen & Kamionkowski 2004). Then, 
N,{z) = [(l + 2/Hc)/(l + 4/Hc)]iVb(2) ~ 0.88 iVb (2), where 
/He = 0.0789 is the helium-to-hydrogen number ratio. 

We assume that the energy transfer through ionization 
results in a photon loss, whereas Compton scatterings reduce 
the average energy of the photons without changing their 
number. So, for z < z' the equations for the cosmological 
evolution of n{z,z') and E(z,z') are 

n{z + dz,z') _ , , p, '^.^dt 

n{z,z') - ^-'l>.^o.{z,Eiz,z))-dz (8) 

E{z + dz,z') g,/ '\\d,t dz . , 

___ = i_^,,_(,,£;(,,,))_d,+ __(9) 

where the energy equation also keeps into account the cos- 
mological redshifting of photons (also note that dt/dz is neg- 
ative). 
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These equations need to be supplemented with the in- 
jection of new photons at each redshift z = z': 



Ci fiomiz') ^dz 

dz 



n{z' + dz, z) 

E{z' + dz,z') = C2mDMC^ ^1 + ^ 



(10) 

(11) 



where nDM(«') is the rate of decrease of the number of DM 
particles per baryon at redshift 2', and mDM is the mass 
of a DM particle (sec Section 2.3). ("i and 1^2 arc numerical 
coefficients which depend on the considered DM particle and 
on the details of its decay or annihilation. 



2.2 Pair production 

Even if the loss of kinetic energy from electrons and 
positrons can be treated exactly in the same way, the anni- 
hilation probability is negligible for the electrons, but must 
be kept into account for the positrons. For this reason, 
when DM decays/annihilations result in the production of 
an electron-positron pair, it is useful to distinguish between 
electrons and positrons. 

In the case of annihilations wo neglect other processes, 
such as the internal Bremsstrahlung, which affect t(z) only 
in a minor way (see Appendix A) . 

Note that in the following we will always include the 
rest energy me(? as part of the energy E of the particle. 



2.2.1 Electrons 

Electrons can transfer their kinetic energy to the IGM 
through collisional ionizations and ionizations by in- 
verse Compton up-scattered cosmic microwave background 
(CMB) photons. Here we neglect the energy loss through 
synchrotron radiation, because the inverse Compton mech- 
anism is more efficient by a factor Ucmb /Us (where [/cmb 
and Ub. are the energy densities of the CMB and of the 
magnetic field, respectively), which is S> 1 unless unrealisti- 
cally strong magnetic fields {B ^ 10~* gauss) are assumed 
to exist in the IGM at 2; > 5. 

It is necessary to remark that photons produced by the 
inverse Compton mechanism are not necessarily absorbed 
by the IGM, and we must distinguish between the fractional 
energy loss rate by electrons 



$c- («,£;) = $e,ion(2,£;) + ^e,com{z , E) , 



(12) 



and the fractional energy loss rate actually absorbed by the 
IGM, 



-iz,E) = (/>e,ion(2,-B) + (pe.coiniz , E) 



(13) 



The ionization losses are completely absorbed by the 
IGM, so that 



(2,£)=$e,ion(2,S)^f^ 



{ZiiNiiiz 

+ ^HaiVHe(2) 



In 



©(7) 



i°( "-S''-- )+"M 

\ th,He / 



(14) 



} 



where v and e are the velocity and charge of the elec- 
tron, 7 = E/{meC^) is the electron Lorentz factor, -Eth.H ~ 
13.59 eV (£th,He = 24.6 eV) is the hydrogen (helium) ioniza- 
tion threshold, Zh = 1 {Zue = 2) is the hydrogen (helium) 



atomic number, Nii{z) = Nb{z)/{1 -|-4/He) is the H number 
density, Nueiz) = Nh{z) fne/i^ + 4/He) is the He number 
density, ©(7) = 4, - (| - ^) ln2 + i (l - i)' , and 



imax.H — —7 



27 rriTTmeV 



TJig +27771 e mil 

27^(4mHj^metJ^ 



[max,He — m^ + (4mH)2+27me(4mH) 

(cfr. Longair 1992, and also Lang 1999, and Chen & 

Kamionkowski 2004). 

The total Compton fractional losses arc given by 



^e,coni{z,E) = ^— ^ gCrTCOradTcMB(2)'' (7^ - 1) 



(16) 



where arad — 7.56 x 10~ ergcm" K~ is the Stefan- 
Boltzmann constant for the radiation energy density, and 
Tcmb{z) ~ 2.726 (1 + z) K is the temperature of the CMB 
radiation at redshift z (cfr. Rybicki & Lightman 1979; Lon- 
gair 1992); the fractional loss rate that is actually absorbed 
by the IGM is 

.,2 



0e,com(2, E) = *e,com(2, E) 



7 



72-1 



1/(7) 



^/.5.(rcMB(.)) 

aradCTcMB(«)^ ^ ^ 



where Bi, (T) is the Planck function for black-body radiation 
at temperature T, and 

3 Eth,H 



z/(7) 



(18) 



4 ftp 7^ 

where hp is the Planck constant. The correction to $e,com 
introduced in equation (17) amounts to neglecting photons 
with pre-intcraction frequencies below ^^(7). This is neces- 
sary because on average the post-interaction energies of such 
photons are below Sion.H, so that they are haxdly absorbed 
by the IGM^ 

As in the case of photons, we assume that ionization en- 
ergy losses lead to the disappearance of electrons, whereas 
Compton losses reduce the average energy of the surviving 
electrons. Therefore, the equations describing the evolution 
of the energetic electrons coming from decays and annihi- 
lations are quite similar to those we used for photons. For 
redshifts z < z' we have 



= l-^<,,ion{z,Ee-{z,z'))^dz (19) 

= l-^e,com{z,E^-{z,z'))^dz, (20) 



ne-{z + dz,z') 

ne-{z,z') 
Ee-{z + dz,z') 

Ee-{Z,Z') 

whereas the injection of new electrons is described by 

nc-{z' + dz, z') = Ci '"dm(z') 

Eg- {z' + dz,z') = C2 "IDM c^- 



(21) 
(22) 



where, again, Ci and ^2 depend on the details of the consid- 
ered decaying or annihilating particle. 

2.2.2 Positrons 

In addition to the energy loss mechanisms of electrons, decay 
produced positrons can also annihilate with thermal elec- 

^ Lyman a opacity generally results only in the scattering of the 
photon, rather than in its absorption, and only a small fraction 
of the energy is absorbed (Purlanetto & Pritchard 2006) . 
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trons in the surrounding gas. The fractional energy loss due 
to the annihilations of positrons of energy E is 

*ann = i; iVe (^;) (Tann (E) (23) 

where the annihilation cross-section Uann is (cfr. Beacom & 
Yiiksel 2006) 

3(7T 



8(7 + 1) 



(24) 



Every annihilation emits two photons, each of energy 
E~/ ~ i(i? + rUeC^), as it involves a positron of energy E 
and an electron whose energy is likely to be close to rrJeC^. 
Such photons are absorbed only if the optical depth they 
encounter is sufficiently high. Here we do not follow their 
radiative transfer in detail, and simply assume that the frac- 
tional energy loss which actually goes in the IGM is 



.n(z, E) = $an„(«, E) /l [1 - e""^ (^■■^■^)] 



with 



T^{z,E^) = 



h 

H{z) 



^{z,E^ 



(25) 



(26) 



where H{z) is the expansion rate of the universe at red- 
shift z, and 4>'y,coin{z,Ej) is the Compton fractional energy 
loss of a photon of energy E-y, as defined in equation (7). 
Equation (26) represents a fraction /2 of the optical depth 
encountered by a photon of energy Ej emitted at rcdshift 
z and traveling an Hubble radius, assuming that the bary- 
onic density does not vary. The parameters /i = 0.91 and 
/2 = 0.6 have been chosen in order to maximize the agree- 
ment between a full radiative transfer treatment and our 
simple approximation. 

Equation (25) implies that annihilation energy loss is 
particularly efficient for positrons with very low kinetic en- 
ergy, as <Tann(-B) oc (7— 1)"^^^ whcu E « TWe (and 7 tends 
to 1). 

For this reason, it is convenient to separate positrons in 
two different groups: "fast" positrons (which effectively lose 
energy through ionizations, inverse Compton, and annihila- 
tions), and "thermal" positrons (whose kinetic energy is so 
low that annihilations are their only energy loss mechanism). 

The treatment of fast positrons is very similar to that 
of electrons: their fractional energy loss rate is 

$e+,f (^, E) = *e,ion(^;, E) + $e,com(.?, E) + *an„(2, E), (27) 

and the fractional energy loss rate actually absorbed by the 

IGM is 

(f>e+,f{z, E) = <^o,ion(«, E) + 0o,com(z, E) + (f)^^^(z, E) , (28) 

where the ionization and Compton terms are exactly the 
same as in the case of electrons. 

Their evolution is described by the equations which are 
simple modifications of those given for electrons, taking an- 
nihilations into account 



nt{z + dz, z') 
nt{z,z') 



= 1 - [$e,ion(«,^f(z,2')) 

-F$ann(2,Sf(«,2;'))]^rf« 



Ei{z + dz,z') (TP, '\\dt, 

E,iz,z') = l-^..o^{z,E,iz,z))-dz, (30) 

where the "f" subscripts refer to "fast" positrons (e.g. , Ei is 
the average energy of the fast positrons) . Their injection rate 
and average energy at injection are identical to those given 
in equations (21) and (22) for electrons: n{{z' + dz,z') = 
ne-{z' +dz,z'), E{{z' +dz,z') = E^-iz' + dz,z'). 

The fractional energy losses (total and absorbed by the 
IGM) of thermal positrons are simply 

$e-H,t(2,-B) =$ann(«,£;) (31) 
0e+,t(z,S) =0ann(z,£;); (32) 

their evolution equations are 

nt{z + dz,z') = nt(2:,2:')[l - $ann(z, -Bt(z))] 

+nt{z,z')'i>,s,n{z,Et{z,z)) (33) 
Etiz) = niec' + ^kBTiGM{z) (34) 

where the "t" subscripts refer to "thermal" positrons, fee 
is the Boltzmann constant, and Tigm(2) is the temperature 
of the IGM at rcdshift z. Wc take the injection rate to be 
simply nt{z',z') = 0, as positrons are naturally "fast" when 
they are created. 

In practice, these equations assume that the only mech- 
anism which leads to the disappearance of a positron is its 
annihilation; instead, ionization energy losses simply turn a 
fast positron into a thermal one. 



2.3 The absorbed energy fraction 

The energy injection rate per baryon resulting from the in- 
tegration of equation (3) can be expressed in the form of the 
fraction /abs of the total energy released by the DM which 
is absorbed by the IGM 



fabsiz) 



e{z) 



(35) 



In the case of decaying DM, tidm (which we defined as the 
decrease rate of the number of DM particles per baryon) is 
given by 



nDM(«) = e 

TDM 



t(0)-t(z) 



?^DM,0 
TDM ' 



(36) 



where tidm.o is the number of DM particles per baryon at 
present, tdm is the lifetime of a DM particle, and t{Q) and 
t{z) are the ages of the universe at present and at redshift z, 
respectively. The leftmost equality is valid when tdm S> t(0), 
which is generally the case. 

Instead, in the case of annihilating DM, tidm is 



1 2 3 

h-omiz) ^ -n-oM.,oNh{Q) {ov){l + z) , 



(37) 



(29) 



where (av) is the thermally averaged annihilation cross- 
section; the 1/2 factor is due to two reasons: first, the DM 
is split in half between particles and anti-particles, and this 
needs to be accounted by introducing a correction factor 1/4. 
However, this must be multiplied by 2, as each annihilation 
involves two DM particles. 
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The definition of the absorbed fraction^, /abs, given in 
equation (35) is such that it can be easily plugged into the 
equations commonly used in studies concerned with decay- 
ing and annihilating DM (e.g. Padmanabhan & Finkbeiner 
2005; Zhang et al. 2006). A second advantage of the /abs no- 
tation is that it is independent of uncertain quantities such 
as tdm and (crw). 



3 APPLICATIONS 

The rate of energy absorption per baryon, e(z), that was 
derived in the previous Section, can be used to improve our 
knowledge of the effects of DM decays/annihilations on cos- 
mic reionization and heating. In fact, most of the previ- 
ous works (Chen &: Kamionkowski 2004; Padmanhaban & 
Finkbeiner 2005; MFP06; Zhang et al. 2006) derived only 
an upper limit to these effects, by assuming that all the en- 
ergy emitted during the decay /annihilation is absorbed, or 
considering the energy absorbed fraction, faha{z) as a free 
(and mostly unknown) parameter. 

Now, we have an estimate of e[z) and /abs (2) which 
accounts for all the important physics involved, and we 
can use it to derive the effective influence of DM de- 
cays/annihilations on the cosmic reionization and heating. 
For this purpose, we ran the public version of the code REC- 
FAST (Seager, Sasselov & Scott 1999, 2000), modified to 
account for the energy injection from DM decays and an- 
nihilations. The adopted procedure is mostly the same as 
described in MFP06. Here, we briefly summarize the most 
important assumptions of this method, pointing out the im- 
provements of our present treatment. 

The IGM is heated, excited and ionized by the energy 
input due to DM decays/annihilations. It is important to 
note that the fraction of the absorbed energy going into each 
one of these components is quite unrelated to how the energy 
was deposited in the IGM in the first place. For example, if 
a keV a photon ionizes an atom, the resulting electron will 
generate a cascade of collisions, and the energy of the photon 
will go not only into ionizations, but also into excitations and 
heating. 

In order to treat this process, we assume that a frac- 
tion (1 — a;)/3 (where x is the ionization fraction) of the 
energy absorbed by the IGM contributes to the ioniza- 
tions (Chen & Kamionkowski 2004), and that a fraction 



ified adding the DM energy injection terms: 



r{x) = c 



(1 



(where C=0.9971, a=0.2663 and 



6=1.3163; ShuU & van Steenberg 1985) goes into heating. 
This definition of !F{x) comes directly from a fit to the 

results of the simulations given by ShuU & van Stoonborg 
(1985), replacing the significantly loss accurate form that is 
used in Chen & Kamionkowski (2004), and in MFP06. 

The evolution equations in RECFAST have been mod- 



We remark that the term "absorbed fraction" might be sliglitly 
misleading, as it is theoretically possible to have /abs > 1 in 
scenarios where t{z) is dominated by the absorption from parti- 
cles in the "background", rather than from the ones which were 
produced recently. This might happen, for example, when riDM 
decreases very rapidly. 



t{z) 1 + 4/He 1 - a;H 



V dz 

. / dTiGM \ 



£'th,H 1 + /He 3 

e{z) 1 -I- 4/hc 1 - a^He 
-Eth.Ho 1 + /hc 3 
2e(z) 1 + 4/hc ^ 

3 fee 1 -I- /ho 

^(xh) + /He ^(SHe 



1 + /h 



(38) 
(39) 

(40) 



where Xh (xHe) is the ionized fraction of hydrogen (he- 
lium) atoms, and £ = [H{z) (1 -|- z)]~^ . These equations are 
slightly different from the ones used e.g. in Padmanabhan 
& Finkbeiner (2005) because in our case e{z) is the energy 
absorption rate per baryon, rather than per hydrogen atom. 

We apply this formalism to two different DM candi- 
dates, i.e. sterile neutrinos and LDM, which are expected 
to have the maximum impact on reionization and heating 
(MFP06). In the case of sterile neutrinos only the decay pro- 
cess is allowed (see Section 3.1). For LDM particles (Section 
3.2) we discuss both the decay and the annihilation process. 

We do not consider heavier DM candidates (such as 
gravitinos or neutralinos) because MFP06 have already 
showed that, even assuming /abs = 1, their contribution to 
reionization and heating is completely negligible, and there 
is no point in extending our formalism in order to account 
for them. 



3.1 Sterile neutrinos 

Sterile neutrinos are one of the most popular warm DM 
candidates (Colornbi, Dodelson & Widrow 1996; Sommer- 
Larsen & Dolgov 2001). Their existence is predicted by the 
standard oscillation theory and required by various exten- 
sions of the Minimal Standard Model, such as the z/MSM 
(Shaposlirukov 2006 and references therein). They are mas- 
sive; so they can decay following different channels. In this 
paper we will consider only the so called radiative decay, i.e. 
the decay of a sterile neutrino into an active neutrino and a 
photon. 

The mass of radiativoly decaying sterile neutrinos can 
be constrained by the absence of any detection of X-ray 
lines consistent with photons due to sterile neutrino decays 
in nearby galaxy clusters (Abaaajian, Fuller & Tucker 2001; 
Abazajian 2006; Abazajian & Koushiappas 2006; Boyarsky 
et al. 2006b). Recently, Watson et al. 2006 applied the same 
method to the X-ray emission from the Andromeda galaxy, 
finding an upper limit 



,c < 2.1 keV 



sin ''26' 
lO-'' 



(41) 



where is the mixing angle. This limit is valid for masses 
24 keV. The most stringent constraints for masses 



mi,s c 



2 < 



rriv sC > 24 keV come from the comparison between the 
unresolved X-ray background and the expected contribu- 
tion from sterile neutrino decays (Mapelli & Ferrara 2005; 
Boyarsky et al. 2006a): 



m^sC 



VO.l 



198 



1.55 X 10-11 



(42) 
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For masses lower than ~ 3.5 keV the main constraints 
arise from the positivity of the lepton number. 

On the other hand, the study of matter power spec- 
trum fluctuations provides a conservative lower limit of 
m,ys c? ^2 kcV (Vicl ct al. 2005), even if more recent es- 
timates significantly increase this lower limit [niv s ^ 14 
keV, Seljak et al. 2006; m^^c^ ^ 10 keV, Viel et al. 2006). 
It is worth noting that these lower limits are independent of 
the mixing angle. 

According to the X-ray observational constraints dis- 
cussed above, the minimum possible lifetime for sterile neu- 
trino radiative decays is (Mapelh & Ferrara 2005; MFP06): 



1.000 



TDM = 2.23 X 10'" s 



10 keV 



if 3.5 < m„« c^/keV < 24, and 



TDM = 9.67 X lO^'' s 



25 keV 



6.6 X IQ-^ 
sin 2261 



1.55 X 10" 
sin226> 



(43) 



(44) 



if m,.acVkeV > 24. 

The current number density of sterile neutrinos iVs,o is 
proportional to the current number density of active neutri- 
nos [Na = 3 Afb(0)/(11 J?), whore rj is the baryon-to-photon 
density] . If we assume also that sterile neutrinos account for 
all the DM, we can write the number of sterile neutrinos per 
baryon as (cfr. Mapelli & Ferrara 2005): 



nsfi = 5.88 X 10^ 



Pcrit 



10- 



g cm 



0.198/ 



(8kev) (6.13x10-1") ' ^^^^ 



X 10" 



^® h"^ g cm ^ is the critical density 



where pcrit — 1 
of the Universe. 

Imposing noM.o ~ ns,o and moM = m^s we have all 
the ingredients needed to calculate nDM(^;) through equation 
(36), for each considered mass. 

Each sterile neutrino decay produces one active neu- 
trino and one photon, each of them with an energy ~ 
|mDMC^. Since the active neutrino does not interact with 
the IGM, we only need to consider the photon. Then, we use 
equations (10) and (11) with Ci = 1 and (2 = 1/2, in order 
to get the injection rate and average energy of photons, and 
proceed to the integration of equation (3). 

The resulting absorbed fraction, fs,bs{z), is shown in 
Fig. 1. In the case of complete and immediate absorption 
of the photon energy (i.e. the case studied by MFP06) 
/abs(-z) would bo 0.5, because half of the sterile neutrino 
mass-energy is taken away by the active neutrino. It is clear 
that the complete absorption approximation is pretty good 
at high redshift, especially for low mass sterile neutrinos 
(m^, 3 c'^ 4 keV); whereas at low redshift it fails by a pos- 
sibly large factor. In Appendix B, we provide analytical fits 
to the /abs(z) curves shown above, for 3 < z < 1000. 

Implementing e{z) in RECFAST (equations 38-40), we 
derive the effective influence of sterile neutrinos on the ion- 
ization fraction and IGM temperature (Fig. 2). The Thom- 
son optical depth, Te, shown in Fig. 2 and in the following 
figures, has been calculated by integrating the well-known 
formula: 



a 
o 

42 











0.100 




















0.010 
















0.001 






10 20 30 50 100 200 500 1000 

z 



Figure 1. Absorbed fraction as a function of redshift for sterile 

neutrinos of masses between 2 and 50 keV (solid lines) . The dotted 
lines, mostly superimposed to the solid ones, show the fitting 
functions listed in Appendix B. The decrease in /abs which can 
be seen for z Si 800 for the 15, 25 and 50 keV curves might be an 
artifact caused by our choice of the redshift where the integration 
is started (zmax = 1100), as is discussed Appendix B. 



az— car x jMb{z), 
dz 



(46) 



where we take 02 = 1000, i.e. the low-redshift boundary of the 
last scattering surface, and z\ = 5, that is approximately 
the lowest redshift at which our fits of the absorbed fraction 
arc valid (mostly because of our underlying assumption of a 
largely neutral IGM). 

In Fig. 2, the solid thin line represents the effect of 25- 
keV neutrinos if we assume complete absorption by the IGM 
(/abs = 0.5 at every z), whereas the solid thick line was cal- 
culated using the derived absorbed fraction /abs. As could 
be expected from Fig. 1, the difference between the com- 
plete and the effective absorption cases (i.e. the thin and 
the thick line) increases as the redshift decreases. For in- 
stance, the ionization fraction in the total absorption case 
(thin line) is higher by a factor ~ 4 at 2; = 20, which be- 
comes a factor ~ 24 at « = 5. The IGM temperature in the 
case of effective absorption t{z) (thick line) is reduced by 
a significant factor (~ 235 at « = 5) with respect to the 
total absorption case. The Thomson optical depth is quite 
negligible, always remaining < 10~^. 

In conclusion, accounting for the effective energy ab- 
sorption significantly reduces the effect of sterile neutrino 
decays on reionization and heating, when compared to the 
case of total absorption (see MFP06). 

The other two thick lines reported in Fig. 2 represent 
the effects of sterile neutrinos of 4 (dotted line) and 15 keV 
(dashed), using the derived absorbed fraction /abs- One can 
be surprised by the fact that 4-keV sterile neutrinos have a 
higher impact on ionization and heating than more massive 
neutrinos. This result comes from two different factors. 

First of all, it depends on the fact that less massive 
sterile neutrinos have higher /abs, especially at low redshift. 
As one can see from the behaviour of Xe and Te, at high 
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Figure 2. Ionized fraction (bottom panel), Thomson optical 
depth (central panel) and IGM temperature (top panel) as a func- 
tion of redshift due to sterile neutrinos. The thick lines are ob- 
tained taking into account the effective absorbed fraction (Fig. 1) 
for sterile neutrinos of masses 4 (thick dotted line), 15 (dashed), 
and 25 keV (solid). The thin solid line shows the contribution 
of sterile neutrinos of mass 25 keV, if we assume an absorbed 
fraction /abs = 0.5. The thin dot-dashed line represents, from 
bottom to top, the relic fraction of free electrons, their contribu- 
tion to Thomson optical depth and the IGM temperature without 
particle decays. In the top panel, the thin dotted line represents 
the CMB temperature. 

redshift 25-keV sterile neutrinos give a stronger contribution 
to ionization than 4-keV sterile neutrinos. It is only at z ^i, 20 
that this tendency is reversed. 

The second reason is not 'physical', but it depends on 
the state of the art of observations. In fact, we adopted 
for each sterile neutrino mass the shortest lifetime con- 
sistent with observations. Present-day observational con- 
straints happen to be much stronger for a 15-keV (Watson 
et al. 2006) than for a 25-keV sterile neutrino (Boyarsky et 
al. 2006a), independently of the intrinsic properties of the 
decay process. 

3.2 Light dark matter 

We define as light dark matter (LDM) particles all the DM 
candidates whose mass is between 1 and 100 MeV. Such 
particles have recently become of interest, because they pro- 
vide a viable explanation for the detected 511-keV excess 
from the Galactic centre (Knodlseder et al. 2005). If they 
are source of the 511-keV excess, then their maximum al- 
lowed mass rriLDM should be 20 MeV, not to overproduce 
detectable gamma rays via internal Bremsstrahlung (Bea- 
com. Bell & Bertone 2004). If we consider also the produc- 
tion of gamma rays for inflight annihilations of the positrons, 
this upper limit might become ^ 3 MeV (Beacom & Yiiksel 
2006). 

In principle, LDM can both decay and annihilate, pro- 



ducing photons, neutrinos and pairs. We will treat both 
LDM decays and annihilations, making the assumption that 
the only decay/annihilation products are pairs. This repre- 
sents quite an upper limit, because neutrinos do not interact 
with the IGM and MeV photons have a low probability to 
be significantly absorbed (see Fig. 3, and also the discussion 
in Chen & Kamionkowski 2004). 

As we did for sterile neutrinos, we assume that LDM 
particles compose the entire DM. So the current number of 
LDM particles per baryon is 

n.BM,o = 4.46 X 10 .p.^^g em-3 j (oigs) 

/ mLDMcn "' / iVb(O) y 

\^ IMeV ) \^2.5 X 10-^ cm-3y ' ^ ' 

where mLDM is the mass of a LDM particle. 
3.2.1 Decays 

The LDM lifetime can be derived by assuming that LDM de- 
cays produce the detected 511-keV emission from the Galac- 
tic centre (Hooper & Wang 2004): 

TDM ~ 4 X 10 s j . (48) 

From equations (47) and (48) we can derive h-DM{z), defined 
in equation (36). 

We assume that LDM decays produce pairs. So, the 
parameters needed in equations (21) and (22) are Ci = 1 
(each decay produces a single electron and a single positron) 
and ("2 = 1/2 (both the electron and the positron receive 
approximately half of the available energy). 

Having defined these values, we then found the rate of 
energy absorption per baryon, e{z), by integrating equation 
(3). The corresponding energy absorption fraction /abs is 
shown in Fig. 3, where the cases of mLDM (? = Z and 10 
MeV are shown (dashed and solid line, respectively). 

In the case of pair production, the assumption of im- 
mediate and complete energy absorption corresponds to 
/abs = 1 at every redshift, because both the electron and 
the positron energy can be absorbed. The effective value 
/abs is always significantly less than 1. 

At high redshift [z ^ 100 — 200), the absorbed fraction 
is relatively high (0.3-0.7, depending on the particle mass), 
and it is dominated by the inverse Compton scattering onto 
CMB photons and by the positron annihilation. 

In fact, positron annihilations contribute to /abs only at 
high redshift, because both the annihilation rate (see equa- 
tion 23) and the probability of absorption of the photons 
they produce (equation 26) scale as positive powers of the 
baryon density. 

Furthermore, CMB photons at high redshift are suffi- 
ciently energetic to be scattered up to the ionization thresh- 
old, because {E^) ~ 30 eV(B/5 MeV)^ [(l+z)/1001] (where 
{E.y) is the average energy of the photon after inverse Comp- 
ton scattering and E is the energy of the electron/positron). 
At lower redshift, the starting energies of CMB photons is 
lower, and the energy boost due to the inverse Compton 
is not sufficient to turn them into ionizing photons, so the 
absorbed fraction drops significantly. 
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Figure 3. Absorbed fraction as a function of redshift for LDM 
particles of mass 3 (dashed line) and 10 MeV (solid line) decaying 
into pairs, and for LDM particles of mass 3 MeV (dotted line) 
decaying into photons. 



At low redshift {z ^ 50), the absorbed fraction /abs 
stabilizes, because collisional ionizations become dominant. 
However, /abs (2 < 50) is always ~ 0.1 for 3-MeV particles, 
or ~ 0.03 - 0.04 for 10-MeV particles. We derived a fit for 
the absorbed fraction, reported in the Appendix B. 

For completeness. Fig. 3 also shows the case of LDM 
decays producing an active neutrino and a photon (dotted 
line) for 3-MeV LDM particles. In this case, at high red- 
shift {z ^ 100), the absorbed fraction for the LDM radiative 
decay is quite similar to that for the decay into a pair (es- 
pecially if mLDM ^ 3 MeV). However, at redshifts z ^ 100 
/abs drops to much lower values for photons than for pairs, 
regardless of ruLDM. 

In Fig. 4 we show the effects of LDM decays on reion- 
ization and heating, both considering the energy absorption 
rate e{z) (thick lines) and the upper limit of complete ab- 
sorption (thin lines, /abs ~ 1). For LDM decays the differ- 
ence between the two cases is important, even if less than for 
25-keV sterile neutrinos. For example, if mLDM c^=10 MeV, 
the ionization fraction (the IGM temperature) at 2 ~ 5 is 
a factor ~ 23 (~ 57) higher in the case of total absorp- 
tion than if we consider our estimate of /abs. The Thomson 
optical depth is reduced by a factor ~ 7, and it is only 
Te ^ 2.1 X 10"^ (instead of Te £ 1.5 x 10"^). 



3.2.2 Annihilations 

We now consider the case of LDM particles of mass 
w-LDM =1, 3 and 10 MeV, assuming that they annihi- 
late and produce electron-positron pairs. Recently Zhang et 
al. (2006) found that, in order to be consistent with the 1- 
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Figure 4. Ionized fraction (bottom panel), Thomson optical 
depth (central panel) and IGM temperature (upper panel) as a 
function of redshift due to LDM decays. The thick lines are ob- 
tained using the effective absorbed fraction (Fig. 3) for decaying 
LDM of masses 3 (thick dashed line) and 10 MeV (solid). The 
thin solid (dashed) line shows the contribution of decaying LDM 
of mass 10 (3) MeV, if we assume an absorbed fraction of 1. The 
thin dotted line in the top panel and the thin dot-dashed line in 
all the panels are the same as in Fig. 2. 



yr WMAP results at the 1 a level, the thermally averaged 
annihilation cross-section for LDM annihilations must be'^ 



(cr«) < 2.2 X 10"^^ cm^ s" 



MeV 



(49) 



The absorbed fraction /abs depends on the redshift. How- 
ever, as a conservative approximation, we take the value 
of /abs in equation (49) to be the maximum value /abs, max 
that we derive with our method (see Fig. 5). In particu- 
lar, /abs, max — 0.5 for 1 and 3-MeV LDM particles, and 
/abs, max — 0.9 for 10-MeV particles; these values (and the 
whole function /abs) are actually independent from the value 
we adopt for (crw). 

We chose to use values of {av) which are close to the 
upper limit given by the above formula, i.e. 4, 12, and 24 
xlO"^^ cm^s"^ for mLDM = 1, 3, and 10 MeV, respec- 
tively. Such values of {av) are quite close to those ({o"d)=0.3, 
2.7 and 30 xlO"^^ cm^ s"^, for the same masses) which have 
been inferred by Ascasibar et al. (2006) in order to reproduce 
the 511-keV excess from the Galactic centre. 

From equations (49) and (47) we then derive the rate of 
change of the number of DM particles per baryon through 
equation (37). In this case, the parameters for the injection 
equations (21) and (22) are ("i = 1/2, ^2 ~ 1, because an 
electron and a positron are produced for every annihilation 
(which obviously involves two annihilating particles). 

In eq. (49) the upper limit actually given by Zhang et al. 2006 
was multiplied by a factor of 2 in order to account for differences 
between the two treatments. 
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Figure 5. Absorbed fraction as a function of redshift for annihi- 
lating LDM of mass 1 (dotted line), 3 (dashed line) and 10 MeV 
(solid line). 



The integration of equation (3) with these parameters 
provides us with the energy absorption rate per baryon, 
e(2), and the corresponding absorbed fraction /abs, which 
is shown in Fig. 5. The behaviour of /abs (2) is quite simi- 
lar to the case of LDM decays. At high redshift {z ^ 100) 
fa.ha{z) is close to the complete and immediate absorption 
value (/abs = 1) and it is dominated by inverse Compton 
scattering and positron annihilations. At low redshift coUi- 
sional ionizations alone contribute to /abs, which suffers a 
large drop between the two regimes. 

We note that the contribution of positron annihilations 
to the absorbed fraction is particularly crucial for low mass 
LDM particles. In fact, for the case mLDM = 1 MeV, in- 
verse Compton scattering is not able to produce ionizing 
photons, even at 2: ~ 1000. For this reason, in absence of 
positron annihilation, the absorbed fraction for 1-MeV LDM 
particles (dotted line in Fig. 5) would depend only on coUi- 
sional ionization, and its plot would essentially be a straight 
line from /abs(5) ~ 0.2 to /abs(lOOO) ~ 0.08. 

On the contrary, for the highest mass we consider, 
m-LDM ~ 10 MeV, the high redshift bump is essentially 
due to inverse Compton energy loss. 

The impact of LDM annihilations on reionization and 
heating (Fig. 6) is quite different from the case of LDM 
decays. In fact, LDM annihilations start to contribute both 
to reionization and heating already at very high redshift 
{z ~ 800); but their role remains negligible at low redshift. 
In particular, the ionization fraction becomes ~ 10"^ and 
the IGM temperature at z ~ 10 is much lower than 10 K. 
This mainly is due to the fact that the annihilation rate 
depends on the square of the baryon density (see equation 
37). 

The Thomson optical depth reported in the central 
panel of Fig. 6 is quite high: Te = 0.08 — 0.10 for all the con- 
sidered LDM particles, i.e. close to the best fit (re = 0.09) 
of the 3-yr WMAP data. Nevertheless, the effects of LDM 
annihilations can be very hardly detected by WMAP. The 
reason is that, differently from the relatively rapid and large 
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Figure 6. Ionized fraction (bottom panel), Thomson optical 
depth (central panel) and IGM temperature (upper panel) as a 
function of redshift due to LDM annihilations. The thick lines 
are obtained taking into account the effective absorbed fraction 
(Fig. 5) for annihilating LDM of masses 1 (thick dotted line), 3 
(dashed) and 10 MeV (solid). The thin solid line shows the con- 
tribution of annihilating LDM of mass 10 MeV, if we assume an 
absorbed fraction=l. The thin dotted line in the top panel and 
the thin dot-dashed line in all the panels are the same as in Fig. 2. 



variation of the electron fraction occurring in standard reion- 
ization scenarios, the a;^ evolution produced by LDM anni- 
hilations tracks very closely, albeit at a slightly higher level, 
and for a long time the relic abundance one. This behavior 
dilutes the effects of these extra electrons, making their im- 
print on the CMB spectrum very tiny, as can be appreciated 
from Fig. 7. 

However, the Thomson optical depth produced by LDM 
annihilations should influence the CMB spectra at quite 
high multipoles, as implied by the results of Zhang et al. 
(2006). In Fig. 7 we show how annihilating LDM parti- 
cles of mass 1-MeV modify the temperature-temperature 
(TT), polarization-polarization (EE) and temperature- 
polarization (TE) CMB spectra. The spectra have been sim- 
ulated by running the public code CMBFAST (Seljak & Zal- 
darriaga 1996; Seljak et al. 2003), modified in order to ac- 
count for the effects of DM decays/annihilations. The main 
effects of LDM annihilations are a certain damping in the 
TT peaks, a sensible variation of the EE spectra for / ^ 100 
and some negligible distortions in the TE spectra. However, 
the simulated spectra agree within 1-a with the 3-yr WMAP 
results. Our plot is obtained assuming that all the cosmo- 
logical parameters have the best fit value indicated by the 
3-yr WMAP data. Leaving the cosmological parameters free 
to change, it should be possible to get an even better agree- 
ment between the WMAP data and the simulated CMB 
spectra derived accounting for LDM annihilations (Zhang 
et al. 2006). 
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Figure 7. Temperature-temperature (top panel), polarization- 
polarization (central panel) and temperature-polarization (bot- 
tom panel) spectra. Thick lines indicate the CMB spectrum de- 
rived assuming Thomson optical depth Te = 0.09 and a sud- 
den reionization model (consistent with the 3-yr WMAP data); 
thin lines indicate the CMB spectrum derived assuming re = 0. 
Dashed (solid) lines indicate the CMB spectrum obtained (with- 
out) taking into account the annihilations of 1-MeV LDM parti- 
cles. Open circles in all the panels indicate the 3-yr WMAP data 
(Hinshaw et al. 2006; Pago et al. 2006; Spergel et al. 2006). 



4 CONCLUSIONS 

We have modelled the absorption rate of the energy released 
in the IGM by DM decays/annihilations producing photons, 
electron-positron pairs and neutrinos. Our model suitably 
describes the energy deposition of a wide class of DM candi- 
dates, such as sterile neutrinos and LDM. Useful fits to the 
absorbed energy fraction as a function of redshift for the 
various particles are given in Appendix B. 

In the case of radiatively decaying sterile neutrinos 
(mass 2-50 keV), at z > 200 a fraction /abs — 0.5 of the 
particle energy is transferred to the IGM, predominantly 
via ionizations; at lower redshifts /abs decreases rapidly to 
values of 0.0005-0.3 depending on the neutrino mass. LDM 
particles can decay or annihilate. In both cases /abs ~ 1 at 
high (> 300) redshift, due to positron annihilation and in- 
verse Compton scattering, and it drops to values around 0.1 
below z = 100. 

Our determination of /abs has a dramatic impact on the 
results of previous studies (which adopted naive assumptions 
for this parameter) concerned with the IGM heating by DM. 

To illustrate this point, we have re-calculated the IGM 
thermal and ionization history induced by either sterile neu- 
trinos or LDM particles, using the previous findings. We find 
that sterile neutrino (LDM) decays are able to increase the 
IGM temperature by 2 = 5 at most up to 4 K (100 K). Both 
these values are 50-200 times lower than the estimates based 
on the assumption of complete energy transfer to the gas. 
In addition, significant departures from the adiabatic tem- 



perature evolution induced by the Hubble expansion occur 
only below z ~ 30, at an epoch when heating and ionization 
by conventional sources (stars or accretion-powered objects) 
are likely to swamp the DM signal. 

LDM annihilations instead produce a very extended 
(5 < z < 800) electron fraction plateau, at a level of 5-10 
times the relic one. The main effect of these extra electrons 
is to extend the cosmic time interval during which the IGM 
kinetic temperature is coupled to the CMB one down to 
2 ~ 100. Although the electron scattering optical depth in 
this case is large (0.08-0.10), its effects on the CMB temper- 
ature/polarization spectra are hardly appreciable. 

The detailed computation of the /abs presented in this 
paper and summarized by the fits given in the Appendix B, 
might be useful for a large number of future applications 
in which the cosmological role of the DM is investigated. 
Among these are the effects of DM decays/annihilations on 
the 21 cm emission (Shchekinov & Vasiliev 2006) and on 
the structure formation history (Shchekinov & Vasiliev 2004; 
Biermann & Kusenko 2006; Ripamonti, Mapelli & Ferrara 
2006). 
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APPENDIX A: INTERNAL BREMSSTRAHLUNG 

In Section 2.2 we have considered the energy losses of electrons and positrons produced by DM decays or annihilations. 
However, because of the so called internal Bremsstrahlung process, it is possible that an annihilation results in the production 
of an electron-positron pair, plus a photon (Beacom et al. 2005). As the photon carries away a fraction of the energy released 
by the annihilation, this mechanism could influence /abs- 

Wc can estimate the fraction of energy carried away by the photons produced by internal Bremsstrahlung by using the 
cross-section reported in Beacom et al. (2005) 



~dE 



a 1 

• ftot 

TT E 



In 



- 1 



(Al) 



where E is the energy of the photon, atot is the total annihilation cross-section in the three-level approximation (Peskin & 
Schroeder 1995), a is the fine structure constant, s = 4 (ttidm c^)^, and s' = 4mDM (? (mDM (? — E). 
Then, the fraction of the annihilation energy which is carried away by the photons is 



dE 



2a 

EdE < — ln(2mDM/m.e) 

TT 



(A2) 



(A3) 



2 atot moM 
where we have exploited the fact that 
dffBr , 4:a ln(2mDM/me) 

■ r-. — ''"tot Fi ■ 

dE TT E 

So, even for the highest mass we consider (tudm = 10 MeV), fsr < 0.017. Such a fraction is small; but, if the internal 
Bremsstrahlung photons were completely absorbed, it is enough to significantly alter the value of /abs at low redshift, at least 
for 10-McV LDM annihilations (sec Fig. 5). 

However, this is not the case: the typical energy of a photon produced by the internal Bremsstrahlung mechanism is 
^ ^moMC^, and equation (26) can be used in order to estimate the optical depth encountered by such a photon for each 
Hubble length it travels. For the 5-MeV photons typically produced by the annihilations of 10-MeV LDM particles this is 
T ~ 2.2 X 10~*(1 -|- z)^^^ , so that the complete absorption scenario is realistic only for z ^ 150, when fsi <^ /abs. Then, it is 
justified to neglect the internal Bremsstrahlung process in estimating /abs- 



APPENDIX B: FITS OF THE ABSORBED ENERGY FRACTION 

It can be useful to derive a fit of the energy absorbed fraction (/abs) for each considered DM particle. We calculated fits with 
errors smaller than 5% for 5 ^ z ^ 1000. In the redshift range 3 ^ z ^ 5 the errors are larger ( ^ 10 — 20%); however, at such 
low redshift our absorbed fractions are likely inaccurate by a larger factor, as our assumption of a mostly neutral IGM breaks 
down. 



Bl Sterile neutrinos 



For the case of sterile neutrino with mass 2 < mu,s<? I keV < 10, a fit was found depending only on the redshift and the 
mass of the particle. 



/abs (2, m,., 3) 



0.5-1-0.032 



8 keV 



This general formula does not hold for higher masses, for which we have found 



(Bl) 



/abs(2,15keV) = 0.99 



/abs(^,25keV) = 0.89 



/abs(2,50keV) = 1.08 



0.5 0.032 



/15keV\ 
V 8 keV y 



0.5 0.032 



'25 keV 
. 8 keV 



0.5 -h 0.032 



/ 50 keV y 
V 8 keV y 
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As can be seen in Fig. 1, these last fits do not take into account the decrease of /abs at very high redshift, which is 
caused by the fact that we start integrating equation (3) at redshift 1100. In fact, if we start the integration at z = 1500, 
the high redshift discrepancies between the fitting formulae above and the actual /abs disappear; however, such a procedure 
is uncertain because the equations we use for describing the energy absorption might not be applicable at an epoch before 
recombination. 

Luckily, the energy injection from sterile neutrino decays is completely negligible at such high redshift, and this uncertainty 
by a factor ^ 2 can be safely ignored. 



B2 Light Dark Matter 

For LDM particles decaying into pairs, we derived two different fits, for 3 and 10-MeV particles, respectively. 



U^iz, 3 MeV) = 0.49 (1 + z)-°-^ exp 



+ 



/abs(«, 10 MeV) 



+ 0.058 (1 + 2)"-^* exp 



0.21(1 + 0)-°-^^ exp 



/ 164 \' 
\l + z) 



0.265 exp 



/ 1220 \ 



5.1 

l + z 



+ 0.222(1 + 0)°-^* exp 



Finally, for annihilating LDM particles, we derived the following fits. 

195 ^ 



U.{z, 1 MeV) 
U,{z, 3 MeV) 



0.32 (1 + z)-"-^^ + 0.55 (1 + z)°-°'^ exp 



_ / 195 y 
\l + z) 



X exp 
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V 1900 J 
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+ 2.9xl0-*(l + 2)expf--^') 
\ 1 + zj 



350 
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/abs(2, 10 MeV) = 0.064 (1 + zy°-^^ + 0.335 (1 + ^;)° " exp 
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